Interacting Phantom Energy 


Zong-Kuan Gud 
Institute of Theoretical Physics, Chinese Academy of Sciences, P.O. Box 2735, Beijing 100080, China 


Yuan-Zhong Zhang 
CCAST (World Lab.), P.O. Box 8730, Beijing 100080 
Institute of Theoretical Physics, Chinese Academy of Sciences, P.O. Box 2735, Beijing 100080, China 
(Dated: August 26, 2017) 


We investigate the role of a suitable interaction between a matter fluid and a phantom field for 
the coincidence problem. There exists a stationary scaling solution which is a stable attractor at 
late times. Furthermore, the cosmic doomsday is avoided in one region of the parameter space. 
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Scalar field plays an important role in modern cosmol- 
ogy. Dark energy can be attributed to the dynamics of a 
scalar field, called quintessence fy, which convincingly 
realize the present accelerated expansion of the universe 
by using late-time attractor solutions in which the scalar 
field mimics the perfect fluid in a wide range of param- 
eters. But regarded as dark energy, quintessence field 
with the state equation parameter w > —1 may be not 
consistent with recent observation |3]. In order to ob- 
tain w < —1, phantom field with reverse sign in its dy- 
namical term may be a simplest implementing and can 
be regarded as one of interesting possibilities describing 
dark energy 4. The physical background for phantom 
type of matter with strongly negative pressure would be 
found in string theory [5]. Phantom field may also arise 
from a bulk viscous stress due to particle production 
or in higher-order theories of gravity M, Brans-Dicke and 
non-minimally coupled scalar field theories |8]. The cos- 
mological models which allow for phantom matter appear 
naturally in the mirage cosmology of the braneworld sce- 
nario [d and in k-essence models ag. In spite of the fact 
that the field theory of phantom fields encounters the 
problem of stability which one could_try to bypass by 
assuming them to be effective fields a, it is neverthe- 
less interesting to study their cosmological implication. 
Recently, there are many relevant studies of phantom 
energy : 


The physical properties of phantom energy are rather 
weird, as they include violation of the dominant-energy 
condition, naive superluminal sound speed and increas- 
ing energy density with time. The latter property ul- 
timately leads to unwanted future singularity called big 
rip. This singularity is characterized by the divergence 
of the scale factor in a finite time in future pa. To avoid 
the cosmic doomsday, specific scalar field models were 
proposed u4. It requires a special class of phantom field 
potentials with local maximum. Moreover, the energy 
density of the phantom field increases with time, while 
the energy density of the matter fluid decreases as the 
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universe evolves. Why the energy densities of dark mat- 
ter and phantom energy are of the same order just at 
the present epoch? In this paper we investigate the role 
of a possible coupling of dark matter and phantom field 
for the coincidence problem. With the help of a suitable 
coupling hd, we find that there exists a stationary 
scaling solution and demonstrate numerically that it is 
a stable attractor at late times. Furthermore, the cos- 
mic doomsday is avoided in one region of the parameter 
space. 

We consider the case in which both the phantom en- 
ergy with constant w < —1 and the cold dark matter 
are present. Then, if the universe ceases to be matter- 
dominated at cosmological time tm, the scale factor can 
be written as 
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at t > tm. It is easy to see that the scale factor blows up 
at t = wtm/(1 + w). This occurs because, even though 
the energy densities in ordinary types of matter are red- 
shifting away, the energy density of phantom energy in- 
creases in an expanding universe. It is possible that both 
components decrease with time if there is a transfer of 
energy from the phantom field to the matter fluid. Then 
the cosmic doomsday may be avoided. 

For a spatially flat FRW universe with the matter fluid 
Pm and the phantom field ¢, Friedmann equation can be 
written as 


2_ K? 

H? =~ (py + Pm), (2) 
where K? = 87Gy is the gravitational coupling and the 
energy density and pressure, pp and Pp, of the homoge- 
neous phantom field ¢ are given by 


bp = -18 + V0), (3) 
P, = ~5#-Vi9), (4) 


respectively, in which V (¢) is the phantom field poten- 
tial. We postulate that the two components, pp and pm, 


interact through the interaction term @ according to 


Pp + 3p Pp = =Q, (5) 
Pm + 34 mpm = Q, (6) 
where 
P 
y% = 1+w= Pp ttp (7) 
Pp 
m + Pm 
Ym = (8) 
Pm 


satisfy yp < 0 and 1 < ym < 2. The interaction term Q 
represents an additional degree of freedom which can be 
specified by the existence of solution with a stationary 
energy density ratio r = pm/pp at late times. Using 
Eqs. (8) and (6) we obtain the evolution equation of the 
ratio r: 


par(2 +2 — 3H, +3). (9) 
Pm Pp 


Obviously, the suitable interaction 


TPp 

Q = 3H (m ~ Yo) E (10) 
guarantees the existence of the stationary solution. Note 
that Q > 0, which implies there is a transfer of energy 
from the phantom field to the matter fluid. 

We assume an interaction characterized by Q = 
3H? (pp + pm) where c? denotes the transfer strength, 
which has already been discussed in Ref. fa. When Yp 
and Ym are assumed to be constants, the two stationary 
solutions to Eq. (10) are 


2 
+ Ym — Yp Ym — Yp 
=a ff ee) a. 
f= (5-1) -1, ay 


which imply 0 < r7 < 1 < rf and rīrł = 1 when 
Ym— p > 4c?. A stability analysis of the stationary solu- 
tion indicates that the matter-dominated scaling solution 
rt is unstable, while the phantom-dominated scaling so- 
lution rọ is stable fg. The energy density ratio r evolves 
from the unstable stationary value rf > 1 to the stable 
stationary solution ry < 1, which is clearly compatible 
with the presently favored observational data Qm =œ% 0.3 
and Qa = 0.7. This may provide us with the dynamics 
of the density ratio that is relevant to the solution of the 
coincidence problem. 

Now consider the case of y,, = 1 (i.e. cold dark mat- 
ter). We will investigate the cosmological evolution with 
the stationary density ratio r. Substituting Eq. (0) into 
Eq. (5), it follows that 


; Yp Hr 
3 Hp, =0 12 
Pp + ET Pp , (12) 
which yields 
Pp X a, (13) 


where v = ptr, If the expansion is matter dominated 


until the time tm, then we can write the scale factor as 


m 


ea (1 ae rz) = (14) 


Using the definitions of r and yp, the Friedmann equation 
vic+r 


ØB) gives 
a x exp (A) : (15) 


In terms of ¢, (I3) becomes 


Pp X exp (-- oe) : (16) 


—(14+ w) 


Since Yp is a constant, that is, potential and kinetic ener- 
gies of the phantom field remain proportional, it follows 
that 
y 3(1 + 
V ($) = Vo exp (n) ; (17) 


—(14+ w) 


It is reassuring to find an exponential potential because 
this type potentials arise very naturally in the models of 
unification, such as Kaluza-Klein theories, supergravity 
theories and string theories. Similar result has been ob- 
tained in Refs. uj. These authors started with an 
exponential potential in which there is a free parame- 
ter and then investigated the parameter range for which 
there exists a stable scaling solution. We have first con- 
structed a solution with expected properties and then de- 
rived the corresponding potential. In the coupled phan- 
tom scenario, it is surprising that the signs of the expo- 
nent in Eq.(0 are different in different regions of the 
parameter space, which are determined by the signs of 
the parameter v. We will find that the increasing po- 
tential corresponds to a climbing-up phantom field while 
the decreasing potential corresponds to a rolling-down 
phantom field. 

The recent observations indicate that approximately 
0.3 and 0.7 of the total energy density of the universe 
attribute to dark matter and dark energy, respectively. 
To solve the coincidence problem of our present universe, 
the stable stationary ratio of the energy densities should 
satisfy r7 < 3/7, which corresponds to c? < 0.21(1 — yp) 
(i.e. the regions I and II in Figure 1). In the region I 
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for —3/7 < yp < 0, since 0 < v < 2/3, the energy density 
of the phantom field decreases in an expanding universe 
and the universe accelerates without a big rip in the fu- 
ture. However, in the region I 


<1- 
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FIG. 1: Regions of the (yp, c”) parameter space with ym = 1. 
In the regions I and II, the scaling solution is a stable attrac- 
tor. In the region III, there exist no stable scaling solutions 
with r3 < 3/7. The universe accelerates for ever in the re- 
gion I, while the universe leads to a big rip in the future in 
the region II. 


for —3/7 < yp < 0 and 


for yp < —3/7, v < 0 indicates that the universe accel- 
erates until a big rip occurs at t = (1 — 1/v)tm. We see 
that the energy transfer from the phantom field to the 
cold dark matter prolongs the lifetime of the universe. 

To study an explicit numerical evolution of the energy 
density ratio r and the phantom energy density (or the 
matter energy density), it is most convenient to rewrite 
the evolution equations (9) and () (or Eq.(@) as a set 
of two first-order differential equations with two indepen- 
dent variables r and p, (or pm) 


1 m — 
r= ser (24r t- 5), (18) 
Tr C 
Pp = -3P (per +e), (19) 
/ 2 g 
Pm = — 3pm (+n Cc = <) , (20) 


where the prime denotes a derivative with respect to the 
logarithm of the scalar factor, N = lna. Let us consider 
two points, (—0.2, 0.2) and (—1.4, 0.2) in the regions I 
and II of the (yp, c?) parameter space in Figure 1, respec- 
tively. We choose different initial conditions and follow 
the evolution. We find that an initial energy density ra- 
tio decreases and there exist an attractor curve, which 
corresponds to the stationary scaling solution. It implies 
that the universe evolves from a dark matter dominance 
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to a phantom dominance and the energy density ratio 
of the two components becomes a constant ultimately. 
Furthermore, we find that the stationary energy density 
ratio increases as the transfer strength c? increases when 
Yp is fixed. In Figures 2 and 3, the energy densities of the 
two components redshift away. From Eqs. (6) and (7, 
we note that the phantom field rolls down the exponent 
potential. In Figures 4 and 5, the energy density of the 
phantom field increases with time, while the energy den- 
sity of the cold dark matter decreases initially and then 
increases ultimately with time. From Eqs. (6) and (1%), 
we note that the phantom field climbs up the exponent 
potential. 

We have considered above that the present accelerated 
expansion of our flat FRW universe is driven by an in- 
teracting mixture of a matter fluid and a phantom field 
with w < —1. In the absence of interaction, as shown 
in Ref. Ef there exist no scaling solutions because the 
phantom energy increases while the matter energy de- 
creases with time. With the help of a suitable coupling, 
there exists a stable, stationary scaling solution, which 
requires a transfer of energy from the phantom field to 
the matter fluid. Furthermore, we demonstrate numeri- 
cally that an interaction between cold dark matter and 
phantom field can drive the transition from a matter 
dominance to a phase of accelerated expansion with a 
stationary ratio of the energy densities of the two com- 
ponents. This interacting phantom approach indicates 
a phenomenological solution of the coincidence problem. 
The different regions in the (yp, c?) parameter space lead 
to different fates of the universe. In the region I, the 
phantom field rolls down the exponent potential and the 
universe accelerates without the cosmic doomsday. How- 
ever, in the region II, the phantom energy increases as 
the phantom field climbs up the potential, which lead 
to the divergence of scale factor in the future. Coupling 
of the quintessence field to dark matter may be worri- 
some because of quantum corrections to the quintessence 
potential ice It is valuable to study the stability of 
the phantom potential under quantum fluctuations in the 
case with a coupling to dark matter. 

Note added. While this letter is under review, 
Ref. presents a similar idea. 
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FIG. 2: r-pp with Ym = 1, yp = —0.2 and œ = 0.2 
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FIG. 3: r-pm with Ym = 1, yp = —0.2 and œ? = 0.2 
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FIG. 4: r-pp with ym = 1, Yp = —1.4 and c? = 0.2 
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FIG. 5: r-pm with Ym = 1, yp = —1.4 and œ = 0.2 
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